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The  propagation  of  very  low  frequency  (VLF)  and  low  frequency  (LF) 
radio  waves  in  the  earth-ionosphere  waveguide  is  affected  profoundly  by  the 
earth's  curvature.  In  particular,  for  frequencies  greater  than  about  20  kHz, 
some  modes  are  possible  for  which  the  energy  is  concentrated  near  the  base 
of  the  ionosphere  and  the  field  strength  near  the  ground  is  small.  It  is 
useful  to  think  of  such  modes  as  being  composed  of  waves  repeatedly  reflectec 
at  the  inside  surface  of  the  ionosphere.  Using  analogy  with  sound  waves,  we 
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can  call  these  modes  whispering-gallery  modes.  Since  they  depend  only  very 
slightly  on  ground  conductivity,  these  modes  are  also  referred  to  as  earth- 
detached  modes. 

The  Navy  currently  is  pursuing  efforts  leading  to  a  feasibility 
demonstration  of  a  balloon  gateway  communications  system  that  would  employ 
a  balloon-to-bal loon  cross  link  based  on  the  whispering-gallery  propagation 
modes.  This  report  discusses  some  recent  whispering-gallery  theoretical 
advances  and  simplifications. 
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VLF/LF/MF  WHISPERING -GALLERY  PROPAGATION  STUDIES 

INTRODUCTION 


Around  the  turn  of  the  century,  Lord  Rayleigh*"'5  explained  the  efficient 
propagation  of  sound  waves  around  the  inside  dome  of  St.  Paul's  Cathedral,  in 
London,  giving  rise  to  a  whispering-gallery  effect.  Between  a  speaker  and  a 
listener,  both  located  close  to  the  wall,  an  infinite  family  of  sound  ray 
paths  exist,  each  having  one  or  more  reflections  (almost  at  grazing  incidence) 
off  the  surface  of  the  dome.  Since  the  travel  time  of  the  sound  along  each  of 
these  ray  paths  is  similar,  constructive  interference  and,  hence,  amplification 
takes  place.  In  essence,  an  audible  sound  will  travel  around  the  inside  of  the 
dome  with  exceptionally  low  attenuation. 

Analogs  of  the  whispering  gallery  also  occur  within  the  earth,  where  the 
interface  between  two  layers  guides  seismic  rays  in  the  same  fashion  as  the 
dome  guides  acoustic  rays.^  Only  the  positions  of  the  source  and  receiver 
differ  from  those  in  Lord  Rayleigh's  whispering  gallery,  both  being  located 
outside  the  layer  interface  at  the  earth's  surface.  The  degree  to  which  any 
particular  interface  acts  as  a  whispering  gallery  depends  on  its  ability  to 
guide  significant  amounts  of  energy  along  its  inner  surface. 

The  propagation  of  very  low  frequency  (VLF)  and  low  frequency  (LF)  radio 
waves  to  considerable  distances  is  conveniently  treated  by  regarding  the  space 
between  the  earth  and  the  ionosphere  as  a  waveguide.  Several  authors  have 
found  that  the  least  attenuated  modes  are  affected  profoundly  by  the  earth's 
curvature.  In  particular,  for  frequencies  greater  than  about  20  kHz,  some 
modes  are  possible  for  which  the  energy  is  concentrated  in  a  region  near  the 
base  of  the  ionosphere  and  the  field  strength  near  the  ground  is  small.®-7  It 
is  useful  to  think  of  such  modes  as  being  composed  of  waves  repeatedly  reflec¬ 
ted  at  the  inside  spherical  surface  of  the  ionosphere,  the  rays  being  chords 
of  this  sphere  (figure  1).  Using  analogy  with  sound  waves,  we  can  call  these 
modes  whispering-gallery  modes.  Since  they  depend  only  very  slightly  on 
ground  conductivity,  these  modes  are  also  referred  to  as  earth-detached  modes. 

The  Navy  currently  is  pursuing  efforts  leading  to  a  feasibility  demon¬ 
stration  of  a  balloon  gateway  communications  system  that  would  employ  a  bal- 
loon-to-balloon  cross  link  based  on  the  whispering-gallery  propagation  modes. 

The  basic  phenomenology  of  VLF  and  LF  propagation  is  well  known  but  the 
effects  of  high-altitude  signal  injection  and  reception  are  less  well  defined. 
In  particular,  the  only  experimental  data  available® are  from  a  one-point 
two-frequency  experiment  which  is  insufficient  for  verification  of  the  low 
frequency/medium  frequency  (LF/MF)  whispering-gallery  propagation  mode.  It  is 
the  purpose  of  this  report  to  discuss  some  recent  whispering-gallery  theoret¬ 
ical  advances  and  simplifications.  Throughout  this  report,  we  will  consider 
only  propagation  between  vertical  electric  antennas  over  a  sea  water  path. 
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Figure  1.  Khispering-Gallery  Mode. Description 


PROPAGATION  IN  THE  EARTH- IONOSPHERE  WAVEGUIDE 


The  energy  of  propagating  radio  waves  is  confined  principally  to  the 
shell  between  the  earth  and  the  ionosphere  and  this  space  is  frequently  denoted 
as  the  terrestrial  waveguide.  For  long  waves,  the  effective  waveguide  height, 
h,  is  comparable  to  a  free-space  wavelength,  X,  and  the  characteristics  of  wave 
propagation  are  determined  jointly  by  the  properties  of  the  two  guide  bound¬ 
aries. 

There  are  a  number  of  propagating  modes  with  distinct  cutoff  frequencies 
similar  to  those  present  in  the  microwave  range.  But,  unlike  the  highly  con¬ 
ducting  guides  of  the  microwave  range ,  the  upper  boundary  of  the  terrestrial 
waveguide  is  diffuse  and  a  poor  conductor;  the  finite  conductivity  of  the 
earth's  surface  is  also  important.  In  the  extremely  low  frequency  (ELF) 
range,  h  is  less  than  X  and  only  one  waveguide  mode  propagates.  For  VLF,  h 
exceeds  X  and  there  are  several  propagating  modes.  In  the  LF  range,  the  num¬ 
ber  of  significant  propagating  modes  may  exceed  20. 

Several  field  representations  can  be  used  to  characterize  the  terrestrial 
propagation  of  radio  waves.  The  fields  in  a  uniform  spherical  shell  between 
the  earth  and  ionosphere  can  be  expressed  as  a  summation  of  spherical  harmon¬ 
ics,  which  involves  Legendre  polynomials  and  spherical  Bessel  functions  of 
integer  order  n.  This  series  converges  very  slowly.  The  number  of  terms 
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required  is  of  the  oroer  of  10  ka,  where  k  =  2tt/X  and  a  is  the  earth's  radius. 
Although  this  series  is  directly  applicable  in  the  ELF  range,  its  exceedingly 
slow  convergence  may  preclude  its  use  at  LF. 

The  Watson  transformation  changes  the  series  of  spherical  harmonics  into 
a  residue  or  mode  series,  where  the  fields  are  expressed  using  Legendre  func¬ 
tions  and  spherical  Bessel  functions  of  complex  order  v.  Each  term  of  the 
mode  series  can  be  identified  as  an  azimuthal  wave  propagating  in  the  6  (=  D/a) 
direction  with  a  distinct  phase  velocity  and  attenuation  rate. 

Basically  different  is  the  geometric-optical  series,  where  the  individual 
terms  are  identified  as  a  ground  wave  that  propagates  along  the  surface  of  the 
earth  and  sky  waves  that  reach  the  receiver  after  a  number  of  reflections  from 
the  ionosphere.  This  series  provides  results  identical  to  those  of  the  mode 
series  for  distances  up  to  about  2  Mm.  10  The  wave-hop  propagation  theory^ 
extends  the  geometric-optical  series  to  great  distances  (i.e.,  deep  into  the 
shadow  region).  The  theory  lends  itself  readily  to  physical  interpretation, 
particularly  in  the  propagation  of  LF  pulses.  Wave-hop  techniques  classically 
have  been  used  in  the  LF  band,  principally  for  ground-to-ground  transmissions. 
LF  wave-hop  calculations  for  elevated  antennas*”  recently  have  been  compared*7 
successfully  with  mode-theory  results. 

The  individual  terms  of  the  geometric-optical  series  corresponds  to  a 
sequence  of  waves  that  reach  the  receiver  over  different  propagation  paths. 

For  pulse  transmission,  the  individual  delays  of  the  received  pulses  corre¬ 
spond  to  the  path  lengths  of  the  ground  wave  and  the  various  sky  waves.  How¬ 
ever,  in  the  residue  series,  a  number  of  modes  have  phase  velocities  (v) 
greater  than  the  speed  of  light  (c)  and  there  is  no  one-to-one  correspondence 
between  the  delay  of  the  various  pulses  received  and  the  phase  velocities  or 
the  delay  of  the  individual  modes. 

For  the  upper  VLF  and  LF  ranges  (when  the  transmitter  and  receiver  are 
located  at  or  near  the  earth's  surface),  the  lower  order  modes  (where  v  <  c) 
are  excited  very  weakly.  On  the  other  hand,  the  higher  order  modes  (v  >  c) 
are  highly  attenuated.  Intermediate  modes  with  v  -  c  dominate  in  the  mode  sum 
and  the  order  of  the  dominant  mode  increases  with  increasing  frequency.*® 

Recent  calculations*7  for  100  to  300  kHz  daytime  ground-to-ground  propagation 
have  shown  relatively  little  modal  structure,  indicating  that  only  a  few  modes 
are  required. 

When  the  transmitter  and  receiver  are  elevated  to  considerable  heights 
(30  to  50  km),  the  height  gain  factors  (HGF's)  for  the  lower  order  modes  (v  < 
c)  are  substantial,  even  though  the  excitation  factor  is  very  weak.  Thus,  the 
modified  HGF  (i.e.,  the  product  of  the  excitation  factor  times  the  height  gain 
of  the  transmitter  times  the  height  gain  of  the  receiver)  can  be  comparable  to 
or  greater  than  unity  for  these  lower  order  least  attenuated  modes.  However, 
for  the  upper  LF  frequencies,  the  daytime  propagation  mode  structure  is  con¬ 
siderably  more  complicated  for  the  elevated  transmitter  and  receiver  cases 
than  it  is  for  ground-to-ground  propagation. ^  On  the  basis  of  mode  theory, 
this  would  be  anticipated  since  whispering-gallery  type  modes  would  become 
more  influential  with  terminal  elevation.  On  the  basis  of  wave-hop  theory, 
the  added  structure  would  be  anticipated  since  more  multipath  possibilities 
exist  when  the  terminals  are  elevated.  For  ground-to-ground  transmission , 
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only  one  path  applies  to  a  single  ionospheric  reflection.  However,  when  the 
terminals  are  elevated,  there  exist  four  paths,  or  hops,  linking  transmitter 
and  receiver  which  correspond  to  a  single  ionospheric  bounce. ^ 


GENERAL  THEORY 


"7  19 

The  phase  integral  form  of  the  modal  equation  can  be  expressed  as '  ’  » 

h 

R^  exp  -i2k ^C2  +  2z/a^1/,2dz  +  in/2  =  exp(-i2nm)  ,  (1) 

20 


where 

R.  is  the  Fresnel  reflection  coefficient  for  vertically  polarized  waves 
incident  on  the  ionospheric  boundary  at  an  angle  of  incidence  <t>  ( see 
figure  1, 

k  =  2mA, 

C  =  cos  <t> , 
n 

a  =  earth's  radius,  and 

z  =  -aC2/2. 

0  n 

The  solution  of  equation  (1)  for  the  case  where  the  earth  and  ionosphere 
are  perfect  reflectors  is7 


C2  =  [ 

n  l 


f3it  (n  -  1/4) 


]  -  £ 


The  attenuation  rate,  an,  and  the  phase  velocity,  v^,  of  mode  n  are  given 

a  =  kImC2/2  (5) 

n  n 


v  /C  -  1  =  ReC2/2  .  (4) 

n  n 

There  is  a  simple  ray  interpretation7  of  the  modal  equation,  as  illus¬ 
trated  in  figure  1.  The  integral  with  limits  zQ  and  h  is  equal  to  the  (com¬ 
plex)  electrical  length  of  a  "to-and-fro"  path  between  the  caustic  surface  at 
a  +  zQ  and  the  ionospheric  reflecting  boundary  at  a  +  h.  The  quantity  m/2  can 
be  interpreted  as  the  phase  advance  imparted  as  the  "ray"  grazes  the  caustic 
surface . 
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When  the  caustic  surface  is  sufficiently  elevated  above  the  earth's  sur¬ 
face,  we  would  have  a  pure  whispering-gallery  mode  of  the  type  considered  by 
Lord  Rayleigh  many  years  ago.  Figure  2  is  a  plot  of  the  height,  zg,  versus 
frequency  for  the  first  five  modes  (for  the  case  where  the  earth  and  iono¬ 
sphere  are  perfect  reflectors) .  The  assumed  ionospheric  reflection  height  is 
bO  km.  For  the  first  order  mode,  zg  is  approximately  equal  to  10,  25,  and  40 
km  at  25,  50,  and  100  kHz,  respectively. 

We  also  notice  from  this  figure  that  for  frequencies  of  less  than  approx¬ 
imately  50  kHz  there  is,  at  most,  one  mode  of  the  pure  whispering-gallery  type 
(for  h  =  60  km).  However,  as  the  frequency  is  increased  there  are  more  and 
more  whispering-gallery  modes.  For  example,  at  100  kHz,  three  modes  are  pres¬ 
ent.  Moreover,  for  frequencies  in  the  MHz  range,  there  would  be  so  many  whis¬ 
pering-gallery  modes  that  no  single  mode  would  predominate.  It  is,  then, 
preferable  to  abandon  the  waveguide  mode  treatment  and  to  use  a  ray-tracing 
method. 

In  the  region  between  the  earth's  surface  and  the  height  zQ,  the  field  is 
evanescent.  Consequently,  excitation  of  these  modes  by  a  ground-based  trans¬ 
mitter  is  very  weak.  This  fact  is  borne  out  by  the  behavior  of  the  first  order 
mode  excitation  factor  at  the  higher  VLF  frequencies  and  ionospheric  reflection 
heights,  where  the  phase  velocity  is  less  than  the  speed  of  light  and  the  exci¬ 
tation  factor  decreases  exponentially;  that  is, 


Figure  2.  zg  Versus  Frequency  for  h  =  60  km 
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'  0  /  4  3  / 

\  -  :ti/2(v  .  t  ,i/2  exH-  iv)  •  (5) 

n  "o  n 

where 

tn  =  (ka/:)2/3(-C2) 

and 

yQ  =  (ka/2)2/3(2h/a)  . 

The  first  order  mode  excitation  factor-1  is  plotted  versus  frequency  in 
figure  3  for  the  case  where  the  earth  and  ionosphere  are  perfect  reflectors. 
For  example,  at  a  frequency  of  30  kHz  and  an  ionospheric  reflection  height  of 
80  km,  the  first  order  mode  excitation  factor  is  approximately  -24  dB,  while 
for  h  =  90  km,  Aj  '  -36  dB. 

On  the  other  hand,  while  the  excitation  factor  decreases  exponentially, 
the  HGF's  for  the  transmitter,  Gn(z),  and  receiver,  Gn(z),  increase  exponen¬ 
tially;  that  is, 

Gn(z)  ~  V*  tn/l+Ai  ^n  ‘  z)  exP  (+  J  ln2)  (6) 

and 

Gn(z}  ~  V*  tiAAi(tn  •  £)  6Xp(+  Jln2)  ’  (7) 


where 


z  =  (ka/2) 2/3 (2hj/ a) , 
z  =  (ka/2)2/3(2hR/a) , 
hj  is  the  transmitting  antenna  height, 
hR  is  the  receiving  antenna  height,  and 

A.  (t  -  z)  is  an  Airy  function,  as  defined  by  Miller--1  or  Abramowitz 
1  n  and  Stegun.-® 

For  real  arguments,  this  function  peaks  (=  0.535)  at  -1,  oscillates  for  large 
negative  arguments,  and  becomes  exponentially  small  for  large  positive  argu¬ 
ments. 

Gn(z)  is  plotted  in  figure  4  for  the  case  where  f  =  100  kHz,  h  =  70  km, 
and  the  ionospheric  conductivity  is  10"6  S/m.1®  For  this  particular  situ¬ 
ation,  the  attenuation  rates  for  the  first  five  modes  are  all  equal  to  about 
3  dB/Mm,  while  the  excitation  factors  are  substantially  different;  that  is, 

![x  -  4.3  x  10"8,  A2  ~  7.4  x  10’5,  Aj  -  1.1  x  1CT2,  and  A$  -  1.6.  The  peaks 
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FREQUENCY  (kHz) 


Figure  3.  VLF  First  Order  Mode  Excitation  Factor  Versus  Frequency 
for  the  Case  Where  Both  the  Earth  and  Ionosphere 
Are  Perfect  Reflectors 
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Figure  4.  Typical  HGF's  of  LF  Waveguide  Modes 
(f  =  100  kHz,  h  =  70  km,  a.  =  10~6  S/m) 

and  valleys  of  the  individual  mode  curves  correspond  to  the  maxima  and  minima 
of  the  Airy  function  A^(tn  -  z) .  From  this  figure,  we  see  that  the  height 
gains  for  the  lower  order  modes  can  be  substantial. 

The  modified  HGF  £i.e.,  the  product  of  the  excitation  factor  times  the 
height  gain  of  the  transmitter  times  the  height  gain  of  the  receiver)  can  be 
expressed  as 


HGF  =  A  G  (z)G  (z)  - 
n  n  n 


- 7— A. (t  -  z)A. (t  -  z) 

2(yo  -  tn)>/2  1  " 


(s; 


If  both  the  transmitting  and  receiving  antennas  are  sufficiently  ele¬ 
vated,  the  modified  HGF  becomes  comparable  to  or  greater  than  unity.  For 
example,  referring  to  figure  4,  we  see  that  for  h^  =  hR  =  40  km,  the  height 

gain  for  the  second  order  mode  G2(z)  =  G2(z)  -  150.  Since  A2  -  7,4  x  10-5, 
HGF  -  1.7.  Thus,  while  it  is  extremely  difficult  to  excite  these  earth-de¬ 
tached  modes  from  a  ground-based  transmitter,  the  advantages  of  exciting  the 
least  attenuated  modes  by  elevated  sources  can  be  significant. 


The  vertical  electric  field  strength  produced  by  a  vertical  electric 
antenna  located  in  the  earth-ionosphere  waveguide  can  be  expressed  as18"^1 
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n 

(9) 


where  all  lengths  are  in  meters  and  the  radiated  power,  P  ,  is  in  kilowatts. 

K 

Once  a  particular  ionospheric  conductivity  model  is  assumed,  the  indi¬ 
vidual  parameters  of  equation  (9)  (excitation  factor,  height  gain,  phase 
velocity,  and  attenuation  rate)  can  be  calculated  for  each  mode  and  summed. 
Unfortunately,  in  the  LF  range,  the  number  of  modes  to  be  summed  may  exceed 
20. 


The  Navy  recently  has  funded  the  General  Electric  Company  (GE)  to  estab¬ 
lish  an  adequate  computer  model  to  facilitate  high-altitude  signal  injection 
and  reception  analysis  at  frequencies  from  approximately  15  to  500  kHz  through 
modifications  of  the  existing  Defense  Nuclear  Agency  (DNA)  sponsored  weapons 
effects  on  D-region  communications  (WEDCOM)  computer  model, 24 -30  WEDCOM  has 
both  LF  ray  theory  and  VLF  waveguide  computational  codes  that  can  be  selected 
for  utilization  depending  on  the  particular  frequency  to  be  analyzed.  Detailed 
information  on  the  WEDCOM  models  and  the  extensive  modifications  implemented 
to  facilitate  high-altitude  injection  and  reception  propagation  computations 
can  be  obtained  from  the  GE  TEMPO  personnel  at  Santa  Barbara,  CA.  The  main 
result  of  their  analysis  is  that  LF  communications  over  sea  water  between 
high-altitude  transmitters  and  receivers  is  feasible  in  both  normal  and 
nuclear-disturbed  environments. 


SIMPLIFIED  APPROXIMATE  THEORY 


One  drawback  of  computer  simulations  is  that  they  often  provide  little 
physical  insight  into  the  various  propagation  parameters.  However,  by  employ¬ 
ing  a  combination  of  ray  and  mode  theory,  the  following  simplified  whispering- 
gallery  propagation  equation  can  be  obtained: 


Ey  -  -12.2  +  10  log  P^  +  G  -  5  log  aw 
-  10  log(a  sin  D/a)  -  aD  +  HGF  dBV/m 


(10) 


This  approximate  expression  for  the  vertical  electric  field  produced  by 
a  vertical  electric  antenna  (which  is  a  slight  variation  of  equation  (26)  of 
reference  27)  provides  some  physical  insight  into  the  various  propagation 
parameters . 

The  quantities  that  make  up  this  expression  are  the  radiated  power  PR  (in 
kW) ,  antenna  gain  G,  duct  width  w,  spreading  loss  a  sin  D/a,  dominant -mode 
attenuation  rate  a,  and  modified  HGF.  The  spreading  loss  is  easily  determined 
and  the  maximum  duct  width  is  equal  to  the  distance  between  the  earth  and 
ionosphere  (approximately  60  to  90  km) . 

For  PR  =  1  kW,  G=  1.8  dB,  w=h=60km,  5  log  aw  =  -27.9  dB,  and 


9 


TR  o503 


F  ■  -63.3  -  10  log  (a  sin  D/a)  -  aD  +  HGF  dBV/m/kW  ,  (11) 

or 

E  -51.7-10  log(a  sin  D/a)  -  aD  +  HGF  dBuV/m/kW  ,  (12) 

where  a  and  D  are  both  in  Mm  and  a  is  in  dB/Mm. 

If  w  =  h  =  90  km,  the  constants  in  equations  (11)  and  (12)  become  -69.2 
and  50. S  dB,  respectively.  A  few  comparisons  of  equation  (12)  with  mode- 
theory  computer  calculations  are  presented  in  the  appendix. 

The  most  important  propagation  parameter  to  be  determined  is  the  dominant  - 
mode  attenuation  rate.  Figure  5  is  a  plot  of  the  180  kHz  measured  field 
strength  versus  range. 31,3*.  Fr0m  this  figure,  we  see  that  there  is  a  substan¬ 
tial  difference  between  daytime  and  nighttime  attenuation  rates  at  the  upper 
LF  frequencies. 

The  propagation  of  VLF  and  LF  radio  waves  to  substantial  distances  is 
made  possible  by  the  high  reflectivity  of  the  lower  ionosphere  at  oblique 
incidence.  The  latter  is  due  to  the  relatively  sharp  gradient  of  the  electron 
density  in  the  D-region  of  the  ionosphere.  In  fact,  for  many  purposes,  the 
assumption  of  an  abrupt  lower  edge  of  the  ionized  region  has  permitted  an  ana¬ 
lytical  approach  to  the  problem  that  has  produced  useful  results.  In  the  main, 
these  are  confirmed  experimentally.  However,  a  number  of  discrepancies  have 


Figure  5.  Measured  Field  Strength  for  180  kHz 
Ground-to-Ground  Transmission 
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been  observed  that  suggest  that  the  sharply  bounded  model  is  not  entirely 
adequate . 


During  the  past  20  years,  evidence  from  a  number  of  independent  experi¬ 
mental  approaches  has  indicated  that  the  profile  of  the  electron  density  (and 
ionospheric  conductivity)  in  the  lower  ionosphere  can  be  described  approxi¬ 
mately  by  an  exponential  function  of  height. 


For  horizontal  polarization,  the  reflection  loss  for  an  exponentially 
varying  ionospheric  conductivity  profile  is  a  function  of  frequency,  angle  of 
incidence,  and  ionospheric  conductivity  gradient  i.-l  For  vertical  polari¬ 
zation,  the  relationship  is  more  complex.  Z  is  the  length  over  which  the  con¬ 
ductivity  changes  by  a  ratio  of  2.71;  it  is  that  reciprocal  of  Wait's  S 
factor. -1  (That  is,  the  commonly  employed  daytime  3  value  of  0.3  km-1  corre¬ 
sponds  to  an  Z  of  3.33  km.)  The  magnitude  of  the  reflection  loss  for  vertical 
polarization  is  presented  in  figure  6  as  a  function  of  frequency;^  cos  = 
0.16,  which  is  the  case  for  grazing  incidence  at  the  earth's  surface.  The 
straight-line  portions  of  the  theoretical  curves  are  also  valid  for  horizontal 
polarization.  Typical  values  of  Z  in  the  VLF/LF  range  are  0.5  to  2  km  at 
night  and  2  to  3.33  km  during  the  day. 


It  should  be  noted  that  the  ionospheric  conductivity  gradient  of  impor¬ 
tance  is  the  value  of  2.  at  the  ionospheric  reflection  height  HR,  Typical 
values  of  HR  in  the  VLF/LF  range  are  55  to  70  km  during  normal  daytime  propa¬ 
gation  conditions  and  75  to  90  km  during  normal  nighttime  propagation  condi¬ 
tions.  The  major  reflection  height  usually  can  be  determined  by  examining  the 
detailed  vertical  field  structure2^  ,34-40  (i.e.,  HR  equals  the  height  where 
the  downcoming  horizontal  magnetic  field  strength  component  is  reduced  by  a 
factor  of  2).  A  simpler  method,4!  which  is  actually  quite  accurate  for  many 
cases,  is  that  the  reflection  height  is  located  at  the  height  where  the  ratio 
of  the  ionospheric  conductivity  parameter21  wr  to  w  is  equal  to  yjl  cos2  <p; 
that  is, 


HR  -  H 


Z  Zn 


l  ^2"  U)  COS2  0  J 


(13) 


^2  ui  cos2  <j> 

where  H  is  the  exponential  ionospheric  conductivity  profile  reference  height. 


The  formula  for  the  ionospheric  reflection  loss  for  horizontal  polari¬ 
zation  is21 


| R (  -  exp [- (2tt2£/X)  cos  .  (14) 

As  we  previously  mentioned,  the  reflection  loss  is  a  function  of  frequency, 
angle  of  incidence,  and  ionospheric  conductivity  gradient.  Since  the  whisper¬ 
ing-gallery  mode  attenuation  rate  is  equal  to  the  reflection  loss  (which 
varies  as  cos  <J>)  divided  by  the  distance  between  reflections,  d1  (which  also 
varies  as  cos  <j>)  ,  the  attenuation  rate  is  independent  of  the  angle  of  inci¬ 
dence  and  is  just  a  function  of  frequency,  ionospheric  conductivity  gradient, 
and  radius  of  the  earth.  That  is, 


=  M  = 

2dj  2a  cos  ij> 


(15) 
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Figure  6.  Ionospheric  Reflection  Loss  Versus  Frequency  Ceos  <j> 
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or,  in  dB/Mm, 

-8.686  x  106  f2ir2i\ 

a  -  ^  r~ y — J  ~  1 3 . 4  ( £/ A )  (16) 

Based  on  an  analysis  of  about  30  years  of  10  to  100  kHz  propagation  data, 
Fierce-1-  devised  the  following  empirical  formula  for  the  attenuation  rate  a 
that  had  a  varying  as  the  frequency  to  some  power  u: 


a 


8.686  x  106. 

o  V.^  •  1 


fkHz)U  dB/Mm  , 


(17) 


where,  for  propagation  over  a  sea  water  path,  u  -  0.5  to  0.6  at  night  and  0.9 
to  1.0  during  the  day. 


We  will  now  take  the  attenuation-rate  formula  that  was  obtained  from  the 
reflection  loss  (equation  (16))  and  modify  it  so  that  both  the  frequency  and 
gradient  vary  as  some  power  u;  that  is,  let 


a 


13.4(Jt/X)U 


dB/Mm  . 


(18) 


If  we  take  typical  daytime  values  of  l  =  2.5  km  and  u  =  0.9,  the  resulting 
daytime  attenuation-rate  formula 


ctD  -  O.lSfj^9  dB/Mn  (I9) 

is  almost  identical  to  Pierce's  empirical  result^  of 

a  -  0.17f°*9  dB/Mm  .  (20) 

^  kHz 

If  we  take  a  typical  nighttime  value  of  l  =  1  km,  the  resulting  nighttime 
attenuation-rate  formula  is  also  very  close  to  Pierce's  empirical  result. 42 

We  will  now  compare  our  attenuation-rate  formula  (equation  (18))  with  the 
results  obtained  from  some  GE  computer  simulations , ^ » ^7 ,28  por  the  compari¬ 
son  presented  in  figure  7,  the  slope  of  the  daytime  ionospheric  conductivity 
profile  was  2.5  km  while  the  nighttime  gradient  was  0.5  km.  The  circles  are 
the  computer  simulation  results  while  the  curves  are  the  attenuation  rates 
determined  from  equation  (18).  Following  Pierce, 42  u  was  assumed  to  be  0.9 
during  the  day  and  0.6  at  night.  Note  that  there  is  excellent  agreement 
between  the  approximate  dominant-mode  attenuation-rate  formula  and  the  com¬ 
puter  simulation  results. 

Figure  8  is  a  similar  comparison  for  the  case  where  the  daytime  gradient 
employed  was  3.33  km  while  the  nighttime  gradient  was  1  km.  Again,  the  agree¬ 
ment  is  excellent.  Note  that,  as  the  frequency  is  increased,  there  is  a  sub¬ 
stantial  difference  between  daytime  and  nighttime  attenuation  rates.  For 
example,  at  400  kHz  the  daytime  attenuation  rate  is  approximately  50  dB/Mm 
while  the  nighttime  attenuation  rate  is  only  15  dB/Mm. 
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Figure  7.  Comparison  of  Attenuation  Rate  Versus  Frequency 
(Daytime  l  =  2.5  km  and  Nighttime  l  =  0.5  km) 


Figure  8.  Comparison  of  Attenuation  Rate  Versus  Frequency 
(Daytime  i  =  3.33  km  and  Nighttime  l  =  1.0  km) 
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A  comparison  of  the  modified  HGF  versus  frequency  is  presented  in  figure 
9.  (The  transmitting  and  receiving  antennas  are  assumed  to  be  located  at 
heights  around  30  to  40  km.)  These  approximate  values  were  determined  with 
the  help  of  the  computer  simulations.  Here,  again,  we  see  substantial  differ¬ 
ences  between  daytime  and  nighttime  propagation  at  the  higher  frequencies. 

For  example,  at  400  kHz,  the  HGF  is  only  about  10  dB  at  night  but  is  near  30 
dB  during  the  day.  However,  this  daytime  HGF  is  easily  offset  by  the  substan¬ 
tial  daytime  attenuation  rate  of  40  to  SO  dB/Mm.  Moreover,  since  many  whis¬ 
pering-gallery  propagation  modes  will  be  present  at  the  higher  LF  and  lower  MF 
frequencies,  the  model  interference  pattern  will  be  severe. 

Since  we  have  now  determined  the  various  inputs  to  the  simplified  whis¬ 
pering-gallery  propagation  equation,  we  will  put  it  to  a  test.  An  experiment 
was  conducted” >9  in  May  1972  that  consisted  of  low  power  (approximately  50  W) 
transmissions  on  219  and  438  kHz  between  balloons  spaced  1.6  Mm  apart  hovering 
for  over  30  hr  at  altitudes  of  30  to  37  km.  Ground-to-ground  as  well  as  bal- 
loon-to-bal loon  signals  were  simultaneously  measured  for  comparison  of  path 
loss.  The  balloon-to-balloon  measured  field  strengths  were  40  to  60  dB  higher 
than  the  ground-to-ground  measured  field  strengths. 

Figure  10  shows  a  comparison  of  the  measured  and  predicted  balloon-to- 
balloon  field  strengths.  The  spread  in  the  predicted  value  is  due  solely  to 
the  variation  in  the  ionospheric  conductivity  gradient  assumed.  The  higher 
daytime  predictions  correspond  to  a  gradient  of  2  km  while  the  lower  predic¬ 
tions  are  for  l  =  3.33  km.  The  higher  nighttime  field  strength  predictions 
are  based  on  a  gradient  of  0.5  km  while  the  lower  predicted  value  assumes  l  = 
2.0  km.  Since  the  measured  and  predicted  intervals  overlap,  we  see  that  the 
measured  field  strengths  can  be  explained  satisfactorily  in  terms  of  realistic 


Figure  9.  Approximate  HGF  Versus  Frequency 
(h^,  =  h^  =  30  to  40  km) 
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Figure  10.  Comparison  of  Measured  and  Predicted  Field  Strengths  for 
the  Balloon-to-Balloon  May  1972  Whispering-Gallery 
Propagation  Experiment 


ionospheric  conductivity  parameters.  However,  the  large  spread  also  points 
out  the  need  to  measure  simultaneously  the  ionospheric  conductivity  profile 
during  a  whispering-gallery  propagation  experiment. 

A  plot  of  the  predicted  nighttime  field  strength  versus  range  for  the 
case  where  the  gradient  is  1  km  is  presented  in  figure  11.  The  attenuation 
rate  varies  from  3  dB/Mm  at  25  kHz  to  16  dB/Mm  at  400  kHz.  The  radiated  power 
is  1  kW  and  autumn  nighttime  noise  conditions4^  are  assumed.  The  circles  rep¬ 
resent  the  range  where  the  signal-to-noise  ratio  (SNR)  is  0  dB  in  a  1  Hz  band¬ 
width.  Note  that  this  range  varies  from  10.5  Mm  at  25  kHz  to  5  Mm  at  400  kHz. 

Figure  12  is  a  plot  of  the  predicted  daytime  field  strength  versus  range 
for  the  case  where  the  ionospheric  conductivity  gradient  is  2.5  km.  The 
attenuation  rate  varies  from  3.5  dB/Mm  at  25  kHz  to  40  dB/Mm  at  400  kHz. 

Again,  the  radiated  power  is  1  kW  and  autumn  daytime  noise  conditions4-^  are 
assumed.  Note  that  the  0  dB  SNR  (in  a  1  Hz  bandwidth)  range  varies  from  11  Mm 
at  25  kHz  to  only  3  Mm  at  400  kHz.  These  last  two  figures  clearly  point  out 
the  advantages  of  employing  frequencies  of  less  than  100  kHz  for  long-distance 
whispering-gallery  propagation. 
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Figure  11.  Predicted  Nighttime  Field  Strength  Versus  Range 
for  Various  Frequencies  =  1  km) 
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Figure  12.  Predicted  Daytime  Field  Strength  Versus  Range 
for  Various  Frequencies  (£  =  2.5  km) 
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DISCUSSION 


The  upper  V'LF  and  lower  LF  frequencies  (-  25  to  50  kHz)  not  only  have  the 
advantage  of  lower  attenuation  rates  but  also  the  advantage  of  very  few  whis¬ 
pering-gallery  modes,  in  some  cases  only  one.  Thus,  the  whispering-gallery 
propagation  model  interference  pattern  will  be  nowhere  near  as  severe  as  it  is 
at  the  higher  LF  and  lower  MF  frequencies. 


The  principal  unknown  in  the  earth-ionosphere  waveguide  model  is  the 
ionospheric  conductivity  which  is,  in  part,  a  function  of  the  electron-  and 
^ion-densitv  distribution  with  height.  Fixed  location  sounding  systems  are 
v  somewhat  inaccurate  because  the  ionization  densities  are  quite  low  in  the 
Region  of  importance  to  VLF  and  LF  propagation. 


Propagation  studies  using  long  path  VLF/LF  data  provide  an  indirect  but 
more  useful  description  of  the  ionosphere  for  propagation  prediction.  Much 
*  success  has  been  achieved  with  such  studies  by  using  an  exponential  electrons- 
onlv  ionospheric  profile  specified  by  a  gradient  i  (=  1/S)  and  a  reference 
height  H.  On  the  basis  of  the  analysis-^  of  the  latest  experimental  data, 
i  =  2.0  km  and  H  =  70  km  during  summer  daytime  propagation  conditions  while 
Z  =  3.53  km  and  H  =  74  km  for  daytime  propagation  in  the  winter. 


Figure  13  is  a  plot  of  the  predicted  30  kHz  ground-to-ground  propagation 
field  strength  versus  range  for  summer  daytime  propagation  conditions  (l  =  2.0 
km,  H  =  70  km).  The  individual  curves,  labeled  1,  2,  and  3,  are  the  individual 
mode  results  while  the  curve  labeled  1,2,3  is  the  resultant  multiple-mode  field 


Figure  13.  Predicted  30  kHz  Ground-to-Ground  and  55  km-to-55  km 
Field  Strengths  Versus  Range  for  Typical  Summer  Daytime 
Propagation  Conditions  (5.  =  2  km,  H  =  70  km,  PR  =  1  kW) 
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strength.  These  curves  were  calculated*^  employing  the  following  values  for 
the  various  propagation  parameters  (as  given  by  Wait  and  Spies'^6):  cij  =  2.54 
dB/Mm;  a,  =  4 . 07  dB/Nlm;  a3  =  9.S5  dB/Ntm;  A,  =  -10.4  dB;  \2  =  2.6  dB ;  and  A  = 
l.S  dB. 

Since  the  third  order  mode  is  highly  attenuated,  it  is  essentially  negli¬ 
gible  beyond  about  5  Mm.  However,  the  first  and  second  order  modes  are  com¬ 
parable  at  all  ranges,  resulting  in  a  substantial  interference  pattern, 
especially  in  the  S  to  10  Mm  range. 

Referring  to  Wait  (pp.  379-389)  we  see  that  for  f  =  30  kHz,  h  =  70  km, 
and  hy  -  35  km,  the  HGF  for  the  first  order  mode  is  (z)  -  2.05  (6.25  dB) . 

On  the  other  hand,  the  HGF  for  the  second  order  mode,  G2(z),  is  minimal  at  hy 
-  35  km.  If  the  ionosphere  were  a  perfect  reflector,  G2(z)  would  be  zero  at 
hy  ~  35  km.  For  1  =  2  km,  G2(z)  -  0.2  (-14.0  dB)  at  this  height.  Therefore, 
if  both  the  transmitting  and  receiving  antennas  are  located  at  a  height  of 
approximately  55  km,  the  modified  HGF's  are 

(HGF) j  =  A1G1(z)G1(z)  -  -10.4  +  6.25  +  6.25  =  +2.1  dB  (21) 

and 

(HGF)2  =  A2G2(z)G2(z)  -  +2.6  -  14.0  -  14.0  =  -25.4  dB  .  (22) 

Also  presented  in  figure  13  (the  upper  and  lower  curves)  are  the  field 
strength  versus  distance  plots  for  the  first  and  second  modes  for  the  situa¬ 
tion  where  both  the  transmitting  and  receiving  antennas  are  located  at  a 
height  of  approximately  35  km.  These  plots  were  calculated  using  both  the 
mode-theory  and  simplified  approximate-theory  equations  (equations  (9)  and 
(12)).  Both  methods  yielded  nearly  identical  results  (within  1  dB)  at  all 
ranges . 

Since  this  first  order  pure  whispering-gallery  mode  is  clearly  dominant, 
the  model  interference  pattern  will  be  negligible  beyond  approximately  2  Mm. 
The  35  km-to-35  km  field  strength  is  3  to  30  dB  greater  than  the  ground-to- 
ground  field  strength  (depending  on  range) ,  with  the  average  enhancement  being 
approximately  12.5  dB  for  ranges  greater  than  9  Mm. 

It  should  be  noted  that  Bahar1*?  has  recently  computed  the  VLF  electromag¬ 
netic  fields  along  the  propagation  path  in  an  irregular  spheroidal  model  of 
the  earth- ionosphere  waveguide  using  a  full -wave  approach.  He  has  also 
pointed  out  the  advantages  of  exciting  the  VLF  earth-detached  mode  to  minimize 
the  effects  of  mode  interference. 


CONCLUSIONS 


In  this  report  we  have  discussed  some  recent  theoretical  advances  and 
have  developed  a  simplified  approximate  theory  dealing  with  whispering-gallery 
(or  earth-detached  mode)  propagation  in  the  earth-ionosphere  waveguide.  We 
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have  shown  that  the  simplified  approximate  theory  agrees  remarkably  well  with 
mode-theory  computer  calculations . 

We  have  also  shown  that  the  bal loon-to-bal loon  measurement  results 
obtained  during  the  May  1 9 2  whispering-gallery  experiment  can  be  explained  in 
terms  of  realistic  ionospheric  conductivity  profiles. 

We  have  determined  that  substantial  HGF's  can  be  obtained  at  the  higher 
LF  and  lower  MF  frequencies.  However,  they  are  offset  by  very  high  attenua¬ 
tion  rates  and  severe  modal  interference  patterns. 

We  have  shown  that  the  upper  VLF  and  lower  LF  frequencies  (approximately 
25  to  50  kHc)  are  more  favorable  for  whispering-gal lery  propagation.  They  not 
only  have  the  advantage  of  lower  attenuation  rates  but  also  the  advantage  of 
very  few  whispering-gallery  modes,  in  some  cases  only  one.  Thus,  the  whisper¬ 
ing-gallery  modal  interference  pattern  will  be  nowhere  near  as  severe  as  it  is 
at  the  higher  frequencies. 

We  have  also  shown  that  during  typical  summer  daytime  propagation  con¬ 
ditions,  field  strength  enhancements  of  5  to  30  dB  (depending  on  range)  can 
be  obtained  by  properly  exciting  the  whispering-gallery  propagation  mode. 
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Appendix 


COMPARISON  OF  SIMPLIFIED  APPROXIMATE  THEORY  KITH 
MODE -THEORY  COMPUTER  CALCULATIONS 


It  is  the  purpose  of  this  appendix  to  compare  some  field  strength  calcu¬ 
lations,  using 


Ev  '  51.7  -  10  log (a  sin  D/a)  -  aD  +  HGF  dBuV/m/kW  ,  (12) 


with  those  calculated  by  using  computer  programs  based  on  waveguide  mode 
theory. 


Figures  A-l  and  A-2  are  two  26.1  kHz  daytime  field  strength  versus  dis¬ 
tance  curves  from  Moler  and  Bickel.48  The  mode  constants  (eigenangles ,  atten¬ 
uation  rate,  phase  velocity,  and  excitation  factors)  were  calculated  as  a 
function  of  distance  and  direction  of  propagation  from  the  transmitter  using 
the  Naval  Ocean  Systems  Center  (NOSC)  developed  MODESRCH  computer  program. 

(For  background  and  specific  details,  see  references  17  and  49  through  53.) 

The  MODESRCH  program  input  variables  are  the  ionospheric  electron  densities 
and  collision  frequencies  appropriate  to  the  time  of  day  and  location,  the 
surface  conductivity,  the  geomagnetic  dip  angle,  flux  density,  azimuth  with 
respect  to  the  direction  of  propagation,  and  the  radio  frequency. 

Figure  A-l  is  a  plot  of  the  typical  daytime  field  strength  as  a  function 
of  distance  for  transmitter  heights  of  70,000  to  150,000  ft.  Figure  A-2  pre¬ 
sents  the  vertical  electric  field  strength  at  sea  level  as  a  function  of  dis¬ 
tance  for  both  east-west  (EW)  and  west-east  (WE)  propagation  in  a  spread  debris 
nuclear  environment.  The  transmitting  antenna  is  located  at  a  height  of  80,000 
ft.  Because  the  nuclear-perturbed  ionosphere  is  isotropic,  this  curve  repre¬ 
sents  the  field  strength  with  distance  variation  for  all  directions  of  propa¬ 
gation  . 


For  typical  daytime  propagation  conditions,  the  ionospheric  conductivity 
gradient  l  ~  3  km. 33  For  the  more  lossy  disturbed  environments,  a  more  repre¬ 
sentative  value  of  i  is  4  km.  54  Employing  these  values  of  i  in 


a  '  1 3 . 4  ( 2./  X ) U 


dB/Mm 


(18) 


with  u  =  1  results  in  26.1  kHz  attenuation  rates  of  3.5  dB/Mm  for  £  =  3  km  and 
4.7  dB/Mm  for  1-4  km. 


The  circles  in  figures  A-l  and  A-2  are  the  field  strengths  calculated 
from  the  simplified  approximate  theory  (equations  (12)  and  (18)).  Note  that 
the  agreement  between  the  two  sets  of  calculations  is  excellent. 


A-l 


VERTICAL  ELECTRIC  FIELD  STRENGTH  (dByV/m/kW) 


Figure  A-l.  Predicted  Field  Strength  Versus  Range  for  26.1  kHz  Transmissions 
From  70,  80,  90,  110,  and  130  Thousand  Feet 


A-2 


TR  6503 


Figure  A-2.  Predicted  Field  Strength  Versus  Range  for  26.1  kHz 
Transmissions  in  a  Spread  Debris  Nuclear  Environment 
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Figures  A-3  through  A-7  are  Pappert ' s  range  plots*'  for  frequencies  of 
100,  150,  200,  250,  and  500  kHz.  These  calculations  were  made  by  using  the 
lower  extremities  of  the  GE-TEMPO  ambient  day  profile.®®  On  each  plot  are 
curves  for  ground-to-ground  transmission  for  transmitter  and  receiver  at  30  km 
and  for  transmitter  and  receiver  at  50  km.  The  null  in  the  ground-to-ground 
transmission  curves  that,  depending  on  frequency,  falls  in  the  range  from 
about  S00  to  1400  km  is  a  manifestation  of  the  groundwave  and  first  hop  sky- 
wave  interference  null.  Even  up  through  300  kHz  the  ground-to-ground  curves 
show  relatively  little  modal  structure,  indicating  that  only  a  few  modes  are 
required  for  that  configuration.  The  mode  structure  is  considerably  more  com¬ 
plicated  for  the  elevated  transmitter  and  receiver  cases. The  number  of 
modes  utilized  ranged  from  a  dozen  at  100  kHz  to  28  at  300  kHz.  Minimum  modal 
attenuation  rates  ranged  from  about  7.7  dB/Mm  at  100  kHz  to  about  19  dB/Mm  at 
300  kHz . 

Also  shown  in  figures  A-3  through  A-7  are  the  simplified  approximate 
theory  ground-to-ground  transmission  results  (equation  (12)).  For  these  cal¬ 
culations,  only  one  mode  has  been  assumed.  With  the  exception  of  the  ground- 
wave  and  first  hop  skywave  interference  null,  the  agreement  between  the  two 
sets  of  calculations  is  remarkably  good. 

The  219  kHz  1.6  Mm  range  field  strength  value  (normalized  to  a  transmit¬ 
ter  radiated  power  of  1  kW)  measured  during  the  May  1972  whispering-gallery 
experiment®^  aiso  plotted  on  the  200  and  250  kHz  plots  (figures  A-5  and 
A-6) .  Note  that  the  measured  value  is  very  close  to  the  30  km-to-30  km  pre¬ 
dicted  values. 

A  comparison  of  the  simplified  approximate  theory  with  Pappert1 s*7  150 
kHz  50  km-to-50  km  results  for  both  the  GE-TEMPO®5  and  Deek's5®  ambient  day¬ 
time  ionospheric  conductivity  profiles  is  shown  in  figure  A-8.  The  HGF 
assumed  is  10  dB  (figure  9)  and  the  single-mode  attenuation  rates  employed  are 
10  dB/Mm  for  the  Deek's  profile  and  11.5  dB/Mm  for  the  GE-TEMPO  profile.55  On 
the  average,  the  agreement  is  quite  good  considering  that  only  one  mode  was 
assumed  for  the  simplified  approximate-theory  calculations. 
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Figure  A-3.  Predicted  100  kHz  Range  Plots  for  Gli-THMPO  Daytime  Profile 


Figure  A-4.  Predicted  150  kHz  Range 
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A-10 


Figure  A-8.  150  kHz  Comparison  for  Deek's  and  GE-TF.MPO  Profiles  (h  =  h  =  50  km) 
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